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Abstract 

Relativistic and unfactorized calculations for the nuclear transparency extracted 
from exclusive A(e, e'p) reactions for 0.3 < Q 2 < 10 (GeV/c) 2 are presented for the 
target nuclei C, Si, Fe and Pb. For Q 2 > 0.6 (GeV/c) 2 , the transparency results 
are computed within the framework of the recently developed relativistic multiple- 
scattering Glauber approximation (RMSGA). The target-mass and Q 2 dependence 
of the RMSGA predictions are compared with relativistic distorted-wave impulse 
approximation (RDWIA) calculations. Despite the very different model assumptions 
underlying the treatment of the final-state interactions in the RMSGA and RDWIA 
frameworks, they predict comparable nuclear transparencies for kinematic regimes 
where both models are applicable. 
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The transparency of the nuclear medium to the propagation of protons is an 
issue of fundamental importance. The nuclear transparency provides a mea- 
sure of the probability that a proton of a certain energy escapes from the 
nucleus without enduring inelastic final-state interactions (FSI). The nuclear 
transparency is a useful quantity for studying nuclear medium effects, and 
in particular, it is very well suited for investigations of the so-called color 
transparency (CT) phenomenon, which predicts a significant enhancement of 
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the transmission of protons through nuclei [1,2] once QCD mechanisms start 
playing a role. 

During the last decade, several investigations of the nuclear transparency have 
been carried out using the A(e, e'p) reaction in the quasi-elastic (QE) regime. 
In this kinematic regime, the impulse approximation (IA), where a quasifree 
single-nucleon knockout reaction mechanism is assumed, has been proved to 
provide a good description of the reaction dynamics. Thanks to the electro- 
magnetic character of the initial-state interactions in an A(e, e'p) process, the 
entire nuclear volume can be probed. 

Nuclear transparency measurements with the A(e, e'p) reaction are available 
for a range of target nuclei. The first experiments were performed at Bates for 
Q 2 « 0.3 (GeV/c) 2 [3], and at SLAC for 1 < Q 2 < 7 (GeV/c) 2 [4,5]. Recently, 
measurements at the Thomas Jefferson National Accelerator Facility (TJNAF) 
provided precise data for the target nuclei 2 D, 12 C and 56 Fe and Q 2 = 3.3, 6.1 
and 8.1 (GeV/c) 2 [6]. The same facility provided an alternate set of data for 
the target nuclei 12 C, 56 Fe and 197 Au and 0.64 < Q 2 <3.25 (GeV/c) 2 [7,8]. 

The prediction of the nuclear transparency to protons poses a serious challenge 
for models dealing with the A(e, e'p) reaction due to the wide range of proton 
energies which are probed in the present-day experiments. As a matter of fact, 
at present there is no uniform and realistic framework in which the proton- 
nucleus FSI effects can be computed for proton kinetic energies ranging from 
0.3 to several GeV. For kinetic energies up to around 1 GeV, most theoretical 
A(e, e'p) investigations are performed within the context of the Distorted- 
Wave Impulse Approximation (DWIA), where the effect of the scatterings 
on the emerging nucleon is estimated with the aid of proton-nucleus optical 
potentials. The parameterizations of these potentials are usually not available 
for proton kinetic energies T p beyond 1 GeV. Beyond this energy, the Glauber 
model, which is a multiple-scattering extension of the eikonal approximation, 
offers a valid and economical alternative for describing FSI. In a Glauber 
framework, the effects of FSI on the A(e, e'p) observables are computed directly 
from the elementary proton-nucleon scattering data through the introduction 
of a profile function. The Glauber method postulates linear trajectories and 
frozen spectator nucleons, and the lower limit of this treatment to A(e, e'p) 
has not yet been established. 

Numerous predictions for the nuclear transparencies within the context of non- 
relativistic Glauber theory have been reported in literature [9,10,11,12,13,14,15,16,17]. 
These results are typically obtained in a non-relativistic and factorized model 
for dealing with the e + A — > e' + (A — 1) + p reaction dynamics. In this con- 
text, non-relativistic refers to the fact that the calculations use bound-state 
wave functions or nuclear densities from solutions to a Schrodinger equation 
and non-relativistic expressions for the electromagnetic photon-nucleus inter- 
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action Lagrangian. In the context of modeling A(e, e'p) processes, factorization 
refers to the approximation of decoupling the electron-proton from the nuclear 
dynamics part in the calculations. 

In this paper we focus on relativistic and unfactorized descriptions of nuclear 
transparencies extracted from quasi-elastic A(e, e'p) processes. In the past, 
relativistic distorted-wave impulse approximation (RDWIA) A(e, e'p) calcu- 
lations for the nuclear transparency have been presented by Kelly [18] and 
Meucci [19]. Kelly used an effective current operator containing the Dirac 
potentials, two-component bound states and distorted waves obtained as so- 
lutions to relativized Schrodinger equations. Meucci used bound-state wave 
functions from a relativistic mean-field approach, while the effective Pauli re- 
duction was adopted to construct the ejectile's wave function. 

For T p > 1 GeV many partial waves need to be computed, thus a description 
of the FSI mechanisms in terms of phenomenological optical potentials may 
not be the most economical way at higher energies. In addition, the global 
optical potentials we use in this work were fitted to data with a limited range 
in T p below 1 GeV, hence the RDWIA studies of the nuclear transparency do 
not cover the kinematic ranges beyond T p —1 GeV. 

Recently, a relativized and unfactorized version of the Glauber model has 
been proposed [20,21,22]. In this paper, transparencies obtained within this so- 
called relativistic multiple-scattering Glauber approximation (RMSGA) frame- 
work will be compared with those obtained in the RDWIA framework as it 
has been implemented by the Madrid-Sevilla group [23,24,25,26]. The com- 
parison is made in a consistent way. For the transparency results which will 
be presented below, this implies that the two frameworks only differ in the 
way they treat the final-state interactions. All the remaining ingredients are 
kept identical. 

In what follows we will first sketch the basic ingredients entering the RDWIA 
and the RMSGA frameworks, thereby indicating the similarities and differ- 
ences between the two. Then, the energy dependence, expressed in terms of 
the four- momentum transfer Q 2 = q 2 — u 2 , and target-mass dependence of the 
nuclear transparencies obtained in the two approaches will be compared and 
confronted with the data. 

Adopting the IA, the central quantity to be computed in a relativistic approach 
to A(e, e'p) is the current matrix element 



where <\>b and 4>p are relativistic wave functions describing the initial bound 
and final outgoing nucleon respectively. Further, J M is the relativistic one-body 




(1) 



3 



current operator. The bound-state wave function <j>B is a four-spinor with well- 
defined parity and angular momentum quantum numbers (/cj,, fit,), obtained 
within the framework of the relativistic independent-particle shell model. It 
may be formally written as 



<t>B{r) 



\if, 



(Or)) 



(2) 



with <3>^(f2 r ) the usual spin spherical harmonics. The quantum numbers of the 
bound-state 4> B reflect the nature of the single-hole state in which the residual 
nucleus is left. 

Within the RDWIA framework, <pp is a scattering solution to a Dirac-like 
equation, which includes scalar (S) and vector (V) global optical potentials 
obtained by fitting elastic proton — nucleus scattering data. The scattering 
wave function, expressed in terms of a partial-wave expansion in configuration 
space, reads 

Mr) = E e- i5 h\lm\s F \jtiYr{V PF )^ K {r) , (3) 

V ltj F nam 1 



where ^(f) are four-spinors of the same form as in Eq. (2), except for the 
fact that they have complex phase-shifts and radial functions. The outgoing 
nucleon momentum, energy and spin are denoted as pp, Ep and sf respec- 
tively. 

In the RMSGA framework, the scattering wave function takes on the following 
form 



(4) 



where <p PF SF is a relativistic plane wave. The entire effect of the FSI mech- 
anisms is contained in the Dirac-Glauber phase Q(b,z), which is an A-body 
operator and reads 



G(b,z) = I] [l-Jdf'\Ur')\ 2 0(z'-z)r(b-b') 
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where the profile function for pN scattering is defined as 
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The parameters 0"*^, (3 p n and e P N depend on the proton kinetic energy and 
are obtained through interpolation of the data base available from the Particle 
Data Group [27]. We wish to stress that Glauber-based models come in very 
different flavours and that the RMSGA formulation used here borrows a lot of 
ingredients from the RDWIA approach to A(e, e'p) processes, except for the 
way of computing the effect of FSI mechanisms on the hit proton, which is 
through-and-through different. 

In order to make the comparisons between the RDWIA and RMSGA trans- 
parency predictions as meaningful as possible, all the ingredients in the A(e, e'p) 
calculations not related to FSI, as those concerning the implementation of rel- 
ativists dynamics and nuclear recoil effects, are kept identical. In particular, 
both pictures use the relativistic bound-state wave functions from a Hartree 
calculation with the Wl parameterization for the different field strengths [28]. 
Further, all the results presented in this work are obtained within the Coulomb 
gauge using the so-called CC2 current operator [29]. For the description of 
nuclear transparencies, the effect of Coulomb distortions has been recognized 
as negligible [18]. Therefore, no attempt has been made to correct for the 
Coulomb-distortion effect. 

We wish to stress that the RDWIA, as implemented by the Madrid- Sevilla 
group, and RMSGA codes adopt very different numerical techniques to com- 
pute the scattering wave functions and the corresponding matrix elements 
of Eq. (1). The Madrid RDWIA code employs a partial- wave expansion to 
solve the Dirac equation for the ejectile. The cylindrical symmetry of the 
Glauber phase of Eq. (5) prohibits any meaningful use of this technique in 
the RMSGA calculations. Instead, the multi-dimensional integrals are com- 
puted numerically. In the limit of vanishing FSI mechanisms, however, the 
two codes should predict identical results. In the Glauber approach this limit 
is reached by putting the Glauber phase of Eq. (5) equal to unity. In the RD- 
WIA picture, the effect of FSI can be made vanishing by nullifying the optical 
potentials when solving the differential equations for the radial wave functions 
determining the outgoing nucleon wave functions. In this so-called relativistic 
plane-wave impulse approximation (RPWIA) limit, the Madrid code is prac- 
tically exact, as has been checked by comparing the partial-wave expansion 
results with the analytical ones [30]. In this plane-wave limit and assuming 
identical input options, the RMSGA-Gent and RDWIA-Madrid codes pro- 
duce differential cross sections with an agreement to better than 5%. This 
discrepancy can be partially attributed to the numerical evaluation of the 
multi-dimensional integrals in RMSGA. This comparison gives us confidence 
about the consistency of the calculations. 

The nuclear transparency provides a measure of the likelihood that a struck 
nucleon with kinetic energy T p escapes from the nucleus. The nuclear trans- 
parency can be extracted from the measured A(e, e'p) differential cross sections 



5 



d 5 a exp (e, e'p) on the basis of the following ratio 

rp I _ $A 3 p m !/S.E m dE m S exp (p m , E m , pp) 




(7) 



°A !^p m dPm I A E m dE m S PW jA(Pm, E m ) 



Here, S exp is the experimentally determined reduced cross section 
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where K is a kinematical factor and a ep is the off-shell electron-proton cross 
section, which is usually evaluated with the CC1 prescription of de Forest [29]. 
The quantities A 3 p m and AE m specify the phase-space volume in the missing 
momentum and energy and are commonly defined by the cuts \p m \ < 300 
MeV/c and E m < 80 MeV. These kinematic cuts, in combination with the 
requirement that the Bjorken variable x = 2 ^ ~ 1, guarantee that the 
electro-induced proton-emission process is predominantly quasi-elastic. For 
example, the effects of two-body meson-exchange and isobar currents, which 
are neglected within the IA, have been shown to be at the percent level for 
quasi-elastic kinematics [31,32]. 

In the above equation, Spwia is the reduced cross section within the plane- 
wave impulse approximation (PWIA) in the non-relativistic limit. The factor 
ca in the denominator of Eq. (7) has been introduced to correct in a phe- 
nomenological way for short-range mechanisms and is assumed to be moder- 
ately target-mass dependent. It accounts for the fact that short-range corre- 
lations move a fraction of the single-particle strength to higher missing ener- 
gies and momenta and, hence, beyond the ranges covered in the integrations 
/ dp m j dE m of Eq. (7). The values for ca which are adopted to extract the 
transparency from the A(e,e'p) measurements are 0.9 ( 12 C), 0.88 ( 28 Si), 0.82 



Theoretically, the nuclear transparencies are extracted from the computed 
relativistic A(e, e'p) angular cross sections for the individual single-particle 
states, according to 



This expression reflects the one used to determine T exp . Indeed, in our ap- 
proach, we obtain the "theoretical" transparencies by adopting identical ex- 
pressions and cuts as in the experiments. Essentially, we replace the measured 
A(e, e'p) angular cross sections by the computed ones. In addition, the integra- 
tion over the missing energy j^ Em dE m has been substituted by a sum over all 



( 56 Fe) and 0.77 ( 208 Pb). 
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occupied shells (J2 a ) i n t ne ground state of the target nucleus. Indeed, the rel- 
ativists Hartree approximation does predict bound-state eigenfunctions with 
a fixed energy-eigenvalue and zero width. When determining the denominator 
in Eq. (9), in our calculations the PWIA limit is accomplished by nullifying all 
sources of FSI mechanisms and neglecting those contributions introduced by 
the presence of negative-energy components in the relativistic bound nucleon 
wave function [30]. 

Transparencies have been calculated for the nuclei 12 C, 28 Si, 56 Fe and 208 Pb. 
All numerical calculations are performed in planar and constant (q, uj) kine- 
matics. The adopted values for q and uj are the central values of the kinematics 
in the A(e, e'p) transparency experiments reported in Refs. [3,4,7,8]. For each 
shell a, the kinetic energy of the outgoing nucleon is calculated by means of 
the relationship T p = uj + e a , where e a is the energy eigenvalue of the cor- 
responding single-particle state. Due to the internal motion of the confined 
protons, the ejected protons emerge in a cone about the transferred momen- 
tum. The boundaries of the cone are restricted by the requirement that the 
"initial " proton momentum \p m \ < 300 MeV/c. 

In the RDWIA calculations, we have employed the global S — V parameteri- 
zations of Cooper et al [33], which provide the best phenomenological optical 
potentials to date. As the highest kinetic energy in these parametrizations is 1 
GeV, RDWIA transparencies are obtained up to four-momentum transfers of 
Q 2 PS 1.8 (GeV/c) 2 . Due to its use of the eikonal approximation, the validity of 
RMSGA becomes questionable when approaching low values of Q 2 . For this 
reason the RMSGA model is not used for calculating transparencies below 
Q 2 ~ 0.6 (GeV/c) 2 . Hence, the kinematic range 0.6 < Q 2 < 1.8 (GeV/c) 2 will 
be covered in both the RMSGA and the RDWIA frameworks. 

First, we investigate the sensitivity of the computed transparencies to the 
adopted parameterizations for the optical potentials. In Fig. 1 results for 12 C 
and 208 Pb are displayed as a function of Q 2 for different optical-potential 
parameterizations contained in Ref. [33]. For 12 C, both the predicted Q 2 de- 
pendence and the value of the transparency depend on whether A-dependent 
(EDAD1/EDAD2) or A-independent (EDAIC) fits for the potentials are se- 
lected. For 208 Pb, the noted differences between the different types of optical- 
potential sets are less pronounced. Within the class of A-dependent param- 
eterizations, the versions EDAD1 and EDAD2 give rise to comparable nu- 
clear transparencies. In the remainder of the paper, the EDAD1 version will 
be used. There are various arguments to motivate this choice. First, the A- 
independent parameterization is only available for a very limited number of 
nuclei, and extrapolation to other nuclei has been discouraged [33]. Second, all 
energy-dependent A-dependent parameterizations in Ref. [33] produce similar 
transparency predictions. Finally, the relativistic transparency calculations by 
Kelly [18] and Meucci [19] employed the EDAD1 parametrization. Adopting 
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the same choice makes easier the comparison between these predictions and 
ours. 

In Fig. 2, the transparencies predicted by the RMSGA and RDWIA models 
are displayed as a function of Q 2 and compared to the world data. The 197 Au 
data are compared to 208 Pb calculations. The RDWIA approach systemati- 
cally underestimates the data by roughly 5 — 10%. The presented RDWIA 
transparency results for 56 Fe and 208 Pb are in better agreement with the data 
than those reported in [7]. The RDWIA transparencies obtained in Ref. [19], 
on the other hand, are rather comparable to ours for low Q 2 , the differences 
increasing for higher values. 

The RMSGA framework reproduces the 12 C data. With increasing target mass, 
the agreement worsens, mostly for the lower values of Q 2 . The quality of agree- 
ment achieved for 12 C and the systematic underprediction of the transparen- 
cies for heavier nuclei, was also a feature of the non-relativistic and factorized 
Glauber calculations of Pandharipande and Pieper [6,16]. A global feature 
of the RDWIA and RMSGA calculations presented here, is that they tend to 
underestimate the measured transparencies. Note that this is at variance with 
the RDWIA calculations of [19]. 

As can be inferred from Fig. 2, the RMSGA framework predicts less absorp- 
tion than RDWIA for a light nucleus like 12 C. With increasing target mass the 
opposite holds true and when approaching the heaviest target nuclei consid- 
ered here, the Glauber framework predicts 5 to 10 percent more absorption. 
The measured Q 2 dependence is reasonably well reproduced by both rela- 
tivistic calculations. For low Q 2 the models reproduce the trend of decreasing 
transparencies. For Q 2 > 2 (GeV/c) 2 , the RMSGA transparencies are close 
to constant, in line with the measured ones and those predicted in typical 
non-relativistic Glauber models. In fact, the modest energy variation of the 
transparency in the RMSGA model is a reflection of the fact that the total and 
elastic proton- nucleon cross sections remain fairly constant once T p > 1.7 GeV. 

In Ref. [18] large discrepancies were observed between the DWIA A(e, e'p) 
transparencies and the ones from Glauber calculations of Nikolaev [10,11]. 
In contradistinction, Fig. 2 indicates reasonably good agreement between our 
RDWIA and RMSGA model predictions for the medium-heavy nucleus 56 Fe 
and modest variations in opposite directions when moving to a lighter or 
heavier nucleus. In Ref. [18] the noted differences between the transparen- 
cies obtained from DWIA and those from the particular Glauber approach of 
Refs. [10,11], are attributed to the fact that the latter adopts a closure prop- 
erty in deriving the expression for the attenuation factor. We wish to stress 
that this approximation is NOT used in the RMSGA formulation of Glauber 
theory. In computing the effect of FSI mechanisms on the the A(e, e'p) cross 
sections, the sum extending over the occupied states a in Eq. (9) is carried 
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out in a similar fashion in RMSGA and RDWIA. 

Investigating the attenuation for each individual shell in the target nucleus 
allows one to study the radial dependence of the FSI mechanisms. In the 12 C 
case, for example, the lsi/2 has spatial characteristics which are very differ- 
ent from the IP3/2 orbit. The attenuation for the individual states represents 
also a more stringent test of the (non-) similarity of the optical-potential and 
Glauber-based models for describing proton propagation through nuclei. In 
Fig. 3, the RMSGA and RDWIA predictions for the attenuation for the in- 
dividual shells in 12 C are compared. These numbers are computed according 
to the definition of Eq. (9) without performing the sum over the states a. 
Obviously, the optical-potential approach predicts more absorption for both 
shells. As expected, both models predict a stronger attenuation for proton 
emission from a level which has a larger fraction of its density in the nuclear 
interior. Again, the results of Fig. 3 illustrate that the proton-nucleus (RD- 
WIA) picture and the proton-nucleon picture (RMSGA) are not dramatically 
different in their predictions. These findings provide us additional confidence 
that the "low-energy" and "high-energy" regime can be bridged in a rela- 
tively smooth manner. Note further that the observed tendency of increasing 
12 C transparencies at low Q 2 , can almost be entirely attributed to the lsi/ 2 
orbital. 

The A-dependence of the nuclear transparencies at various values of the four- 
momentum transfer is studied in Fig. 4. The RDWIA framework reproduces 
the measured A-dependence, while RMSGA slightly overestimates it. Under 
the assumption that the attenuation effect is proportional to the radius of the 
target nucleus one would naively expect that the A-dependence of the nuclear 
transparency can be parameterized as 

T(Q 2 )=c(Q 2 )A- a ^ , (10) 

with a— 1/3. In the work of Ref. [7] it was shown that the dependence of 
T eX p(Q 2 ) on the mass number could be nicely fitted with c(Q 2 ) = 1 and 
a = 0.17 ± 0.04(Q 2 =0.65), 0.22 ± 0.05(Q 2 =1.3), 0.24 ± 0.04(Q 2 =1.8), 0.25 ± 
0.04(Q 2 =3.3), 0.20 ± 0.02(Q 2 =6.8). To guide the eyes these curves are also 
displayed in Fig. 4. 

In conclusion, we have presented for the first time a relativistic calculation for 
the nuclear transparency for the process e + A ^ e' + (A — 1) + p covering 
the wide range of quasi-elastic kinematics in which experiments have been 
performed. An optical-potential approach has been used up to the highest 
kinetic energy (T p « 1 GeV) for which potentials are readily available. Beyond 
that region we gathered our results within the context of a relativized and 
unfactorized Glauber framework. In a medium-Q 2 range, both models have 
been applied and their predictions compared. Both frameworks accomodate 
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relativistic effects in the bound-state and scattering wave functions, as well as 
in the electromagnetic current operator. Despite the very different assumptions 
underlying the description of FSI effects in an optical-potential and Glauber 
based approach to A(e,e'p), their predictions for the nuclear transparency 
and, in general, the effect of attenuation for different single-particle levels, are 
comparable. 
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Fig. 1. The sensitivity of the computed nuclear transparencies in C and Au 
to the adopted choice for the parameterization of the relativistic optical potentials. 
Results of RDWIA calculations with the EDAD1 (solid curve), EDAD2 (dashed 
curve) and EDAIC/EDAIPb (dot-dashed curve) are shown. Data points are from 
Refs. [3] (open squares), [4,5] (open triangles), and [7,8] (solid triangles) 
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Fig. 2. Nuclear transparencies versus Q 2 for A(e, e'p) reactions in quasi-elastic kine- 
matics. The RMSGA (solid lines) are compared to the RDWIA (dashed lines) re- 
sults. Data are from Refs. [3] (open squares), [4,5] (open triangles), [6] (solid circles) 
and [7,8] (solid triangles). 
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Fig. 3. The Q 2 dependence of the computed nuclear transparency for the two sin- 
gle-particle orbits in 12 C as obtained in the RDWIA and RMSGA approach. 
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Fig. 4. The A-dependence of the nuclear transparency at five values of the 
four-momentum transfer Q 2 . The solid (dashed) curves are RMSGA (RDWIA) 
calculations. The dotted curves represent the A~ a ^ ' parametrization, while the 
dot-dashed curve gives A' 1 / 3 . Data are from [7,8] (solid triangles) and [4,5] (open 
triangles). 
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